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chains 
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We study the quasiparticle band structure ol isolated, infinite (HF)oo and (HCl)oo bent (zigzag) 
chains and examine the effect of the crystal field on the energy levels of the constituent monomers. 
The chains are one of the simplest but realistic models of the corresponding three-dimensional crys- 
taUine solids. To describe the isolated monomers and the chains, we set out from the Hartree-Fock 
approximation, harnessing the advanced Green's function methods local molecular orbital algebraic 
diagrammatic construction (ADC) scheme and local crystal orbital ADC (CO-ADC) in a strict 
second order approximation, ADC(2,2) and CO-ADC(2,2), respectively, to account for electron 
correlations. The configuration space of the periodic correlation calculations is found to converge 
rapidly only requiring nearest-neighbor contributions to be regarded. Although electron correla- 
tions cause a pronounced shift of the quasiparticle band structure of the chains with respect to 
the Hartree-Fock resuh, the bandwidth essentially remains unaltered in contrast to, e.g., covalently 
bound compounds. 

PACS numbers: 71.20.-b, 71.15.Qe, 71.20.Ps, 31.15.Ar 
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I. INTRODUCTION 

Hydrogen fluorideisSii^ and hydrogen chloride^ii are 

representatives of hydrogen-bonded crystals. This type 
of bondingSiSiiSiiiiiS is weak and essentially preserves 
the electronic structure of the constituent HF and 
HCl monomers, respectively. Hydrogen fluoride forms 
strong hydrogen bonds, whereas hydrogen chloride forms 
only weak hydrogen bonds. Therefore, the two com- 
pounds are good candidates to study the quasiparti- 
cle band structure of hydrogen-bonded periodic systems. 
The structures of both crystals are very similar at low 
temperatures; the monomers are arranged in terms of 
weakly interacting parallel zigzag chains. The isolated 
infinite chain is considered to be a simple but real- 
istic one-dimensional model for the three-dimensional 
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The isolated (HF)oo and (HCl)oo chains have been 
studied Carefully!i3,i4ii5,i6,i7,i8,i9,20,2i^22 particularly the 

ground state has been examined, e.g., Refs. ive 
a detailed ah initio analysis of the electron correlation 
contributions and the basis set dependence of the bind- 
ing energy. The Har tree-Fo ck band structure of (HF)oo 
is discussed in Refs. and the one of the 

(HCl)oo chain is investigated in Refs.'l3''l4 15 18. Corre- 
lation corrections for (HF)oo have been obtained for the 
center (F point) and the edge (X point) of the Brillouin 
zone by Liegener and LadikiS. who use many-particle per- 
turbation theory in the M0ller-Plesset partitioningifiS, No 
correlation corrections have been determined for (HCl) 00 
so far. 

To carry out accurate ah initio correlation calcula- 
tions for periodic systems, we set out from the wannier 
program^'^'^^ which directly produces Wannier orbitals 
by solving modified Hartree-Fock equations i22i2iSii2^ To 
account for electron correlations, an advanced Green's 
function method, the crystal orbital algebraic diagram- 



matic construction (CO-ADC) scheme for Wannier or- 
bitals22i2& is used. It has been devised recently on the 
basis of the well-established ADC scheme for atoms, 
molecules, and clustersiSLS^ It is a method that approx- 
imates the Feynman-Dyson perturbation series for the 
self-energy and contains sums of certain proper and im- 
proper diagrams to infinite order. The scheme sets out 
from the equations in terms of Bloch orbitals^SiSSiS and 
expresses them in Wannier representation which allows 
one to specify efficient cutoff criteria for the configu- 
ration space. The well-known robustness of the ADC 
method even facilitates to study strong electron correla- 
tions which typically occur when inner valence electrons 
are treatedj^S*^ Particularly, they have been examined 
in series of oligomers of increasing length which model 
infinite chains. Yet in this work, we focus on 
the quasiparticle (or main) states for which a one to one 
correspondence to Hartree-Fock one-particle states can 
be established. Like the one-particle states, they form a 
so-called quasiparticle band structure. ^^'^^ The CO- ADC 
theorji22iS& has been implemented on top of wannier in 
terms of the CO- ADC programJ^ 

This article is structured as follows: The ADC and 
CO- ADC schemes in local representation are introduced 
in Sec. EI whereas Sec. IIIII describes geometries, basis 
sets, cutoffs, and the computer programs employed. We 
devote Sec. lIVI to the discussion of the energy levels of the 
HF and HCl monomers paying particularly attention to 
the impact of electron correlations and their implications. 
The quasiparticle band structure of the infinite chains is 
investigated in SecE Conclusions are drawn in Sec. I VII 



II. THEORY 

In order to determine the energy levels of atoms, 
molecules, clusters, and crystals, we set out from the 
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Schrodinger equation with the full nonrelativistic Hamil- 
tonian neglecting nuclear motions; the ground state is 
obtained in restricted, closed-shell Hartree-Fock approx- 
imation. '^i^2i^i*^2iiM4ii^ From the orbital energies, the 
ionization potentials (IPs) and electron afhnities (EAs) 
follow immediately by Koopmans' theorem(^S*^£i^ the 
energy levels are simply given by their negative value, 
i.e., -IP and -EA. 

Taking electron correlations into account, the energy 
levels now result from the poles of the one-particle 
Green's function (or particle propagator )22ii2il2i^2i^ in 
energy representation; it is expressed, using atomic units, 
in the spectral or Lehmann representation bj^SiSfli^ 
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and depends on the crystal momentum fc,S242iiiS The 
residue amplitudes are y'^\k) — ( | eg | ^'^+^(fc) ) 



and 4")(fc) = («'^-i(-fc)|cg^|^^), where cj^ de- 
stroys the electron in the Hartree-Fock Bloch or- 
bital ipj:^{f s). Here r and s correspondingly denote spa- 
tial and spin coordinates. The negative values of the pole 
positions of Gpq{k, uj) are the IPs and the EAs, /„(fc) and 
An{k), respectively. The infinitesimal 77 > in Eq. 
is needed to render the Fourier transformation between 
time and energy variables convergent. 

The pole strengths or spectroscopic factors are defined 
bji^ \x^\k)\'^ and \yi"'\k)\'^ for IPs and EAs, corre- 
spondingly. They quantify the relative spectral inten- 
sities measured in (inverse) photoelectron experiments 
assuming the sudden approximation. '^"■^'^ Summing over 
all Bloch orbitals yields the orbital independent pole 
strengths which are given 
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The inequalities < |y^"^(fc)|Mx),"^(fc)|2, P^"^ < 1 are 

satisfied. As P^^ = 1 holds for the nth IP and EA, 

respectively, obtained in Hartree-Fock approximation via 

Koopmans theoremf^SiiiiS formula 10) is a measure for 

the one-particle character of the nth state. 

To evaluate the particle propagator the 

self-energy I]rs(fc,a;) with respect to the residual 

interaction — here the difference between the full Hamilto- 

nian and the Fock operator — is introduced via the Dyson 
equation2«'43,49,50^ 

G{k, oj) = G"(fc, oj) + G°{k, uj) S(fc, uj) G{k, oj) , (3) 

where Gpg{k, uj) denotes the Green's function of the non- 
interacting system. The self-energy naturally decom- 
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poses into an energy independent part, the static self- 
energy, and an energy dependent part, the dynamic self- 
energyf^SiS 



S(fc,cj) = S°°(fc) + M(fc,cj) , 



(4) 



with lim M{k,uj) = 0. The dynamic self-energy also 

LJ — >±00 

breaks down into two contributions: 



M{k,uj) = M+{k,uj)+M-{k,uj) 



(5) 



The static and the dynamic self-energies are related 

sr,(fc) = EE^fe,fc..[C,fc%][2^/EGL(fc',^) 

X (S^„(^', uj) + Muvik', uj)) G,rik', lu) dcj 

(6) 

Here Vj^^ p ^ ^j:^ are antisymmetrized two-electron in- 
tegrals^SilL^ and the contour integration runs along the 
real axis and closes in the upper complex ui planei^^*^^*^ 
Formula lO can be evaluated using the Dyson expansion 
method (DEM).^?i^'?i^^^ Due to the relation ® , it is suffi- 
cient to concentrate on the approximation of the dynamic 
self-energy because the static part can be obtained from 
it afterwards. 

The algebraic diagrammatic construction (ADC) 
scheme is used to express the dynamic self-energy in a 
nondiagonal representation, 

M,%{k,uj) = U^\k) - K^{k) - C±(fc)]-i Uf{k) , 

(7) 

the so-called general algebraic form or ADC form. The 
nth order approximation of the ADC is constructed by 
inserting the perturbation ansatz in terms of the resid- 
ual interaction U^{k) = ^^'> {k) + (^) + . . . and 
C^{k) = C^^^\k) + C^'-^\k) + --- into Eq. and ex- 
panding the resulting formula into a geometric series^ 
Comparing the resulting terms with the terms arising in 
the Feynman-Dyson perturbation expansion of the dy- 
namic self-energy, analytical expressions can be found for 
the matrices showing up in Eq. (jT)). Thereby, frequently 
linear combinations of the expressions of a few diagrams 
have to be formed. The construction works because both 
series have the same analytic structure. 

To evaluate the equations for the dynamic self-energy, 
we switch to a representation in terms of generalized 
Wannier orbitals. Such a local representation is par- 
ticularly beneficial to treat crystals with a large unit 
(,gjjj37j38j42j45 Suitable schemes to express the crystal mo- 
mentum dependent ADC form Q using Wannier orbitals 
are discussed in Refs. 20 26. In this work, we exclu- 
sively use the fully translational symmetry adapted form 
in Eq. (29) of Ref-IH 

When changing to Wannier orbitals, one has to keep in 
mind that the Fock matrix is no longer diagonal. The off- 
diagonal matrix elements are most efficiently accounted 
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for by a diagrammatic expansion in addition to the per- 
turbative evaluation of the two-electron interactioniS2iS& 
The resulting approximation schemes are denoted by CO- 
ADC(m,n). This notation indicates that the contribu- 
tions which (partly) involve onc-clcctron interactions are 
treated in order m and the contributions that are exclu- 
sively given by two-electron interactions are treated in 
order n. In this article, we resort to the lowest order CO- 
ADC(2,2) approximation in Eq. (32) of Ref. 26. Then, 
in Eq. ISJ) of this work, M+ (A;, oj) is associated with two- 



particle-one-hole {2plh) configurations and M^^(fc, w) in- 
volves two-hole-one-particle {2hlp) configurations. 

With the translational symmetry adapted Fock ma- 
trix F(k) and explicit expressions for the static and the 
dynamic self-energy, the energy levels of a crystal can be 
determined for the crystal momentum k; they are given 
by the eigenvalues E{k) of the Hcrmitian band structure 
matrix, 



Bik) = 



u-\k) 





I F{k) + Ts°°{k) U+\k) 

U+{k) K+{k) + C+{k) 
^ U~{k) K-{k) + C-{k) J 
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which are obtained by solving the eigenvalue problem 

B{k)X{k) = X{k)E{k), X\k)X{k) = 1 . (9) 

Here X{k) denotes the eigenvectors of B{k). Its ma- 
trix elements, Xp„(fc), represent the residual ampli- 
tudes y^p^\k) or x^\k) in Eq. (jSJ depending on whether 
they are of TV -I- 1 or of TV — 1 type. 

As crystals are infinite, the convergence properties of 
the lattice sums occurring in the matrix elements of B(k) 
have to be examined. Careful analyses of the sums in 
expressions which ar e similar to ours are carried out in 
Refs. l57ll53l5<ififlfilL Although convergent, these sums 
have to be truncated suitably using a cutoff procedure 
because the configuration space, i.e., the number of 2plh 
and 2hlp configurations considered to form B{k), needs 
to be restricted. It can be achieved by realiz- 
ing that the importance of individual configurations for 
the description of the energy levels of a crystal can be 
measured by the magnitude of the one- and the two- 
electron integrals involved. This criterion allows a very 
fine grained, individual selection of configurations accord- 
ing to a given cutoff threshold. The procedure is termed 
configuration selectionf^^^^ 

Having set up the band structure matrix, its partial 
diagonalization is carried out efficiently using a complex- 
Hermitian band- (or block-)Lanczos algorithmi22iS& The 
computer code is based on the complex-symmetric ver- 
sion of Sommerfeld that, in turn, is an adaption of the 
real-symmetric implementation by Meyer >S Band Lanc- 
zos methods are iterative algorithms that utilize the 
product of a matrix, here B{k), with a few vectors to 
transform the original matrix into a band diagonal form 
which subsequently can be diagonalized efficientlyi^SiSSiSi 
The band-Lanczos algorithm usually converges to the 
spectral envelope already after a few matrix times vector 



operations. This corresponds to only moderately sized 
band diagonal matrices i^SiS 

A quasiparticle point of view is assumed for the outer 
valence regime and the lowest virtual states, i.e., we focus 
on those states which are essentially of one-particle char- 
acter with only a moderate adjustment due to electron 
correlationsi22iS The number of quasiparticle bands is 
the same as the number of Hartree-Fock bands and a clear 
association between both can be made. To identify the 
quasiparticle states of a crystal in the spectrum of B{k), 
we resort to the orbital independent pole strength (0) 

(n) 

extracting only states with Pj_ > 0.7. In analogy to 
the energies of one-particle states, the quasiparticle en- 
ergy levels group into bands with respect to the crystal 
momentum, forming a so-called quasiparticle band struc- 
ture. 

A local orbital based ADC scheme for isolated 
molecules can easily be devised setting out from the CO- 
ADC equations in Wannier representation:^*''^- let the 
unit cell of a primitive cubic crystal lattice with a macro- 
scopic lattice constant be occupied by a given molecule. 
Then only Fock matrix elements and two-electron inte- 
grals within a unit cell are nonncgligible. Consequently, 
all states of the crystal are A^o-fold degenerate and the 
band structure matrix has only to be diagonalized at 
the r point to obtain all distinct energy levels of the 
crystali22i2& These levels can be identified with the en- 
ergy levels of the isolated molecule. By this line of 
argument, it is justified to drop the lattice sums and 
lattice vectors totally in the CO- ADC equations to ob- 
tain analytical expressions for the local molecular orbital 
ADC scheme. The dependence on the crystal momentum 
quantum number vanishes and the resulting matrix B is 
now termed energy level matrix. Similar to the nomen- 
clature for the methods based on crystal orbital ADC, 
the corresponding approximative schemes for molecules 
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FIG. 1: (Color) Structure of infinite (HF)oc and (HCl), 
zigzag chains with X = F, CI, respectively. 



are referred to as ADC(m,n). 



III. COMPUTATIONAL DETAILS 
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Hydrogen fluoride and hydrogen chloride monomers 
are diatomic molecules of Coov symmetry. The experi- 
mental values for the internuclear distance are 0.91680 A 
for HF and 1.27455 A for HCliSi Upon crystallization to 
a three-dimensional solid, the HF and HCI monomers 
arrange in long zigzag chains. The structure of the sin- 
gle infinite (HF)oo and (HCl)oo chain [Fig. ^ is deter- 
mined by three parameters, the H — X distance r, the 
X ■ ■ ■ X distance i?, and the angle a — L (HXH) with 
X = F, CI. The following experimental values for the 
parameters are taken: r = 0.92 A, R — 2.50 A, and 
a = 120° for (HF)oo^'*^'?^ and r = 1.25 A, R = 3.688 A, 
and a — 93.3° for (DCl)ooi^ Structural information for 
HCI crystals is not available. Yet HCI and DCl crystals 
have very similar lattice constants and are considered 
to be isomorphous.S In this study, we use the cc-pVDZ 
basis set^^'^^'^^ to represent H, F, and CI atoms in the 
compounds. 

A self-consistent solution of modified Hartree-Fock 
equations is carried out with the wannier program,^'^-^^ 
which directly yields Wannier orbitals.^"-'^''-^''-^^ A finite 
cluster of unit cells is utilized by wannier serving as 
a support for the Wannier orbitals in the origin cell. 
It consists of the origin cell and up to fourth nearest- 
neighbor cells (nine unit cells altogether) for (HF)oo and 
up to third nearest-neighbor cells (seven unit cells al- 
together) for the (HCl)oo chain. To carry out Hartree- 
Fock calculations for molecules, we let each unit cell of 
a one-dimensional lattice with a constant of 1 ^m com- 
prise a monomer. This ensures that interactions with 
neighboring unit cells, which contain periodic images of 
the molecule, are negligible. The Wannier orbitals in the 
origin cell can thus be identified with the orbitals of the 
isolated monomer. 

As WANNIER determines only occupied Wannier or- 
bitals, translationally related projected atomic orbitals, 
so-called crystal projected atomic orbitalS)22i2& have been 
devised to be used as virtual Wannier functions i21 They 
are based on the projected atomic orbitals introduced 
by Pulaj!^ and Saeb0 and PulajsS^ and implemented by 
Werner and co-workers et aZ.iS*2i The wannier program 
produces pseudocanonical Wannier orbitals, i.e., the oc- 
cupied block of the Fock matrix in the origin cell is diag- 



FIG. 2: (Color) Ionization potentials and electron affinities of 
a HF and a HCI molecule where IP > and EA < 0. The 
solid green lines denote Hartree-Fock values (Koopmans' the- 
orem^Si^Si^), while dashed red lines depict ADC(2,2) results. 



onalized. Therefore, in the case of molecules occupying 
a primitive cubic crystal with a macroscopic lattice con- 
stant, only the virtual block in terms of the projected 
atomic orbitals contains off-diagonal matrix elements. 

The WANNIER program22i24 also performs the transfor- 
mation of the Fock matrix elements and the two-electron 
integrals from a representation in terms of one-particle 
basis functions to the Wannier representation. These 
quantities are the only ones which enter subsequent cor- 
relation calculations by the local molecular orbital ADC 
and the local crystal orbital ADC methods. In this study, 
we apply the lowest order approximation, ADC(2,2) and 
CO-ADC(2,2), to molecules and infinite chains, respec- 
tively. To this end, we remove the explicit spin depen- 
dence of the equations involved, following the arguments 
of von Niessen et al.S^ The resulting equations for the 
two methods are implemented in terms of the CO-ADC 
programpS 

Both ADC(2,2) and CO-ADC(2,2) employ an addi- 
tional perturbative expansion of the off-diagonal terms 
of the Fock matrix in local orbitals with respect to their 
canonical orbital based counterparts. This introduces 
an extra approximation using ADC(2,2) in relation to 
ADC (2) whose accuracy has been critically assessed in 
Ref. 20 for a single HF molecule. The outer valence IPs 
and the smallest EAs are found to be well described by 
the correlation method. Yet it is not able to describe 
the strong electron correlations which are observed in 
the inner valence region (F 2s dominated states) of the 
HF monomer in the same way as the ADC(2) scheme 
that, nevertheless, is very inaccurate in this region, too.^" 
This limitation of the method does not imply a constraint 
for this work because we focus exclusively on the main 
states of the systems. 
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Level No. -e IP2.2/EA2,2 IP3/EA3 



3 


20.30 


18.92 


19.64 


4, 5 


17.11 


14.64 


15.73 


6 


-4.99 


-4.70 


-4.67 



TABLE I: Ionization potentials and electron affinities of a 
HE molecule. They are given in Hartree-Eock approxima- 
tion by the negative value of the orbital energies — e; the 
ADC(2,2) values are denoted by IP2,2/EA2,2. The reference 
values IP3/EA3 are obtained with the ADC(3) schemeiSS. All 
data are given in electronvolts. 

IV. HF AND HCl MONOMERS 

Before we turn to the infinite chains in Sec. we 
investigate the ionization potentials (IPs) and the elec- 
tron affinities (EAs) of isolated diatomic HF and HCl 
monomers. This will allow to better interpret the 
band structures of the chains later on. On the ab- 
scissa of Fig. 13 we give the IPs (> 0) and EAs (< 0) 
of the molecules while the ordinate indicates the pole 
strength Q of each state. The data are obtained once 
directly from the Hartree-Fock orbital energies by ex- 
ploiting Koopmans' theorem^Si^2ii& and the other time 
in ADC(2,2) approximation. The numerical values are 
listed in Tabs. Ill and UTI aside from reference IPs and EAs 
which were obtained with the ADC(3) scheme^^ using 
GAMESS-UK-Zfi for the molecular Hartree-Fock calculation 
in conjunction with the ADC program of TaranteUiJIWe 
do not present data which correspond to the IPs of core 
states of F Is type and inner valence states of F2s type 
for HF. Neither do we show the core states of CI Is, CI 2s, 
and CI 2p types in conjunction with the inner valence 
states of CI 3s origin for HCl. Above all, only the EA 
with the lowest energy in the given basis set is displayed. 
Yet the resonance properties'^ of the molecules are not 
investigated here further. 

Many of the characteristics of the IPs and the EAs of 
the two molecules can already be understood in the inde- 
pendent particle model. The outer valence region of the 
two molecules is formed by two distinct IPs in the range 
of 14-22 eV for HF and of 10-18 eV for HCl. By inspect- 
ing the molecular orbitals, we find the lowest IPs of both 
monomers to correspond to an ionization from the two 
equivalent 7r-type lone pairs on fluorine and chlorine, re- 
spectively. Consequently, they are twofold degenerate. 
Ionization from the third cr-type pair requires more en- 
ergy as it is oriented towards the hydrogen atom and thus 
is attracted by its positive partial charge. 

In the two spectra, the differences between fluorine 
and chlorine atoms basically manifest themselves in two 
effects. Firstly, there is an overall shift of all IPs and 
EAs of HCl to lower energies with respect to correspond- 
ing IPs in HF because in chlorine the nuclear charge 
is shielded additionally by the lower lying CI 2s and 
CI 2p shells which are missing in the fluorine atom.»24 



Level No. -e IP2,2/EA2,2 IP3/EA3 



7 


16.89 


16.32 


16.38 


8, 9 


12.83 


12.21 


12.30 


10 


-4.06 


-3.48 


-3.44 



TABLE II: Ionization potentials and electron affinities of a 
HCl molecule. The symbols are chosen as in Tab. 13 All data 
are given in electronvolts. 



Secondly, the larger internuclear distance in HCl re- 
sults in a weaker interaction between hydrogen and chlo- 
rine atoms compared with the interaction between hy- 
drogen and fluorine atoms and causes, in conjunction 
with the lower electronegativity of the chlorine atom, the 
HCl monomer to be weaker and more covalently bonded 
than the HF monomer. 

As soon as one accounts for electron correlations, the 
IPs and EAs of HF and HCl in Fig. |21 change consider- 
ably; they shift to lower absolute values with respect to 
the corresponding Hartree-Fock values. '^-^ The displace- 
ment of the IPs in the outer valence region in ADC(2,2) 
approximation is noticeably smaller for HCl compared 
with HF; the contrary holds for the lowest EA. Inspecting 
Tabs, in and HTl one observes that the ADC(2,2) method 
already provides in most cases a considerable improve- 
ment over the Hartree-Fock approximation. However, 
note that Koopmans' theorem value corresponding to en- 
ergy level 3 of the HF molecule is slightly closer to the 
reference values than the ADC(2,2) value. This obser- 
vation is not very surprising because the HF molecule 
is well known to be difficult to describe accurately with 
quantum chemical methods. 

The fact that the outer valence IPs of HCl in Hartree- 
Fock approximation are closer to the ADC(2,2) val- 
ues compared with HF leads to the conclusion that 
electron correlations are less important for these states 
of HCl. The reverse trend holds for the lowest EA of 
the molecules. This finding can be ascribed to the fact 
that the chlorine atom is notably bigger than the fluorine 
atom. Hence, the valence electrons of HF are squeezed 
to a smaller volume, leading to a worse performance of 
Hartree-Fock theory; it is caused by an insufficient de- 
scription of the Coulomb hole around the electrons which 
is especially important in the case of large geometrical 
orbital overlaps.^2i^Si^ 

Electron correlations increase the energy splitting of 
the IPs of HF with respect to the Hartree-Fock result, 
whereas the splitting remains nearly the same for HCl. 
Moreover, the approximately uniform shift of the outer 
valence IPs of HF and HCl indicates that electron cor- 
relations for these states are dominated by the fluorine 
and the chlorine atom, respectively, and are essentially 
unaltered upon formation of the molecular bond. 
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FIG. 3: (Color) Band structure of a (HF)oo chain. Hartree- 
Fock bands are given by the solid green lines. The dashed red 
lines depict CO-ADC(2,2) quasiparticle bands which are de- 
termined by accounting for 2plh and 2hlp configurations that 
involve Wannier orbitals in the origin cell and the nearest- and 
the next nearest-neighbor cells. 



V. HF AND HCl CHAINS 

Having understood the nature of the IPs and EAs of 
the HF and HCl molecules in the previous Sec. IIVI their 
changes upon crystallization can now be studied. 

The Hartree-Fock band structure of the (HF)oo and the 
(HCl)oo chain is displayed in Figs. Oland^l respectively; 
a close-up of the valence bands is shown in Figs. and 
Numerical data of the energy levels of both compounds 
at the F and the X point are given in Tabs. rmi and lIVI ~ 
As a unit cell of the chains comprises two monomers, the 
number of Hartree-Fock bands is twice the number of 
Hartree-Fock IPs and EAs of the monomer [Fig. [5] and 
Tabs, m and ITT]. For (HF)oo, the four low lying occupied 
bands, which mainly correspond to F Is core states and 
F 2s inner valence states, are left out. Similarly, we do 
not show the energy bands originating from CI Is, CI 2s, 
and C\2p core states and the inner valence of CI 3s char- 
acter. Furthermore, only the four lower conduction bands 
in the cc-pVDZ basis seliSLSiSS are given for both com- 
pounds. Again an examination of the resonance proper- 
ties^^ is not pursued further. 

The Hartree-Fock bands which are situated energet- 
ically around — 18eV for (HF)oo and around — 12eV 
for (HCl)oo are formed by the two equivalent 7r-type 
lone pairs of mostly F 2p and CI 3p character of the iso- 
lated molecules at —17.11 eV and —12.83 eV, respectively 
[Fig. Hand Tabs, d and lU. Note that -IP or -EA, cor- 
respondingly, gives the energy level. By inspection of 
Fig. 121 we see that the degeneracy of the molecular or- 
bitals is lifted only slightly due to the crystal field upon 
formation of the infinite chain. The lowest valence bands 
cross in Figs. El and |S1 which is allowed by rod group 
symmetry. 

The energy bands of (HF)oo in the range of — 24eV 



FIG. 4: (Color) Band structure of a (HCl)oo chain. The sym- 
bols are chosen as in Fig.|3 



to -20 eV and of (HCl)oo between -18 eV and -17eV 
are constituted by the third outer valence orbital of the 
isolated monomers which is also predominantly of F 2p 
and CI 3p character with an orbital energies of —20.30 eV 
and of — 16.89 eV, respectively [Fig. El and Tabs. jH and 
lll| . The dispersion of these two energy bands is much 
larger than the dispersion of the four bands discussed in 
the previous paragraph. This indicates that hydrogen 
bonding in the chains is mainly mediated by the orbitals 
of the monomers which lead to these two outer valence 
bands. 

The Hartree-Fock band structure of the (HF)oo chain is 
reported by Berski and Latajkaii for a series of basis sets. 
Among these, the 6-31-f G(d,p) basis is of most compara- 
ble quality to the cc-pVDZ basisSLSiSS employed in this 
study; the plot of the band structure of (HF)oo in Fig. 3 
of Ref. agrees very well with the plot in Fig. of this 
article. We read off the energy of the top of the valence 
bands at the F point from Fig. 3 of Berski and Lata- 
jka;^^ it lies at — 18.3eV which is very close to our value 
of — 17.9eV from Tab. Ivl The Hartree-Fock band struc- 
ture of the (HF)oo chain has also been studied by Liegener 
and Ladiki^ who use a double-^ basis set [(9s 5p) / [3s 2p] 
for fluorine and (6s Ip) / [2s Ip] for hydrogen to which 
they refer as DZP basis set. They observe that unit cells 
up to third nearest-neighbors are sufficient to represent 
the Fock matrix. My computations confirm this finding. 
Liegener and LadikiS. obtain an energy of — 17.57 eV for 
the top of the valence bands and 3.57 eV for the bottom of 
the conduction bands, i.e., their Hartree-Fock band gap 
amounts to 21.14 eV. The top of the valence bands is only 
by 0.3 eV smaller than our result at the F point. How- 
ever, the bottom of the conduction bands is by 1.78 eV 
smaller than ours. The deviation of our result from those 
of Liegener and Ladiki^ can most likely be ascribed to the 
different basis sets employed, i.e., cc-pVDZ versus DZP 
because the DZP basis set lacks a d-function on fluorine 
which is present for cc-pVDZiSLSiSS Furthermore, they 
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FIG. 5: (Color) Valence band structure of a (HF)oo chain. 
Zoom into Fig. |21 



use a slightly different geometry for the (HF)oo chain. 

Blumen and MerkellS. treat the (HCl)oo chain using a 
(3s) / [Is] basis on hydrogen and a (10s 6p) / [3s 2p] basis 
on chlorine. They include up to second nearest-neighbor 
cells to represent the Fock matrix. Their band struc- 
ture [Fig. 1 in Ref. agrees nicely with Fig. 0] in this 
work. We read off from the figure of Blumen and Merkel 
the value — 13.4eV for the top of the valence bands and 
6.9 eV for the bottom of the conduction bands at the 
r point, yielding a band gap of 20.3 eV. My corre- 
sponding values are -13.186 eV and 4.035 eV [Tab. EH 
amounting to a band gap of 17.221 eV which is somewhat 
smaller than the value of Blumen and Merkel. Again this 
can most likely be ascribed to the difference in the one- 
particle basis sets. 

In order to investigate the effect of electron correla- 
tions on the band structures of the chains, we utilize 
the CO-ADC(2,2) method. The convergence of several 
key quantities with respect to the number of unit cells 
whose Wannier orbitals are considered to form 2plh and 
2hlp configurations in correlation calculations is summa- 
rized in Tabs. IVl and I VII for (HF)oo and (HCl)oo, respec- 
tively. Namely, the top of the valence bands, the bottom 
of the conduction bands, the band gap, and the width 
of the upper and the lower valence band complexes are 
shown. They are compared with the plain Hartree-Fock 
data (with "0" unit cells for the configuration selection) . 
In first place, the quasiparticle band structure of the 
chains is obtained considering only a minimum fraction, a 
single unit cell, in correlation calculations.^"" This causes 
upward shifts of the top of the valence bands by 1.90 eV 
for (HF)oo and by 0.54 eV for (HCl)oo, whereas the bot- 
tom of the conduction bands shifts downwards by 0.22 eV 
for (HF)oo and by 0.60 eV for (HCl)oo. Subsequently, 
we enlarge the configuration space to include also con- 
figurations which extend to the nearest-neighbor cells. 
This causes additional shifts of the top of the valence 
bands by 0.16 eV for (HF)oo and by 0.10 eV for (HCl)oo. 
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FIG. 6: (Color) Valence band structure of a (HCl)oo chain. 
Zoom into Fig. 0] 



Moreover, shifts of the bottom of the conduction bands 
of 0.12 eV for (HF)^^ and of 0.15 eV for (HCl)oo are ob- 
served. Further inclusion of second and third nearest- 
neighbor cells has only a minute effect. The other quan- 
tities in Tabs . FVl and I VII exhibit a similarly quick conver- 
gence. 

The fully converged quasiparticle band structures of 
the chains — configurations in up to second nearest- 
neighbor cells are considered — are shown aside from the 
Hartree-Fock band structures in Figs.|21to|ni We observe 
pronounced corrections of the valence bands of (IIF)oo 
due to electron correlations. Their upward shift is much 
larger than the downward shift of the conduction bands. 
In (HCl)oo the contrary holds, however; the value and the 
difference of the shifts are much smaller compared with 
those of (IIF)oo. These observations match the behavior 
of the IPs and EAs found in the monomers in Sec. IIVI 
Moreover, an analysis of the ground-state binding ener- 
gies of the chains shows a similar behavior of short-range 
contributions of electron correlations which are also sig- 
nificantly larger in (IIF)oo than in {}iC\)oo^^^ Yet the 
impact of the configurations from more than a single unit 
cell in correlation calculations was found to be bigger 
in (HCl)oo. 

The observed long-range correlation corrections of the 
quasiparticle band structures of the chains from more 
distant unit cells are, firstly, caused by the adjustment 
of the many-particle system to the extra charge of the 
electron which is added to or removed from it in the one- 
particle Green's function. Secondly, corrections due to 
van der Waals interactions take effecti^ 

Valence and conduction bands of both compounds do 
not shift symmetrically, i.e., by the same amount, up- 
wards or downwards [for numerical values see Tabs. Ivl 
and I VII and the discussion of the convergence of the CO- 
ADC(2,2) calculations above]. Instead, in (HF)oo, the 
valence bands show a much stronger infiuence due to elec- 
tron correlations than the conduction bands. The reverse 
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TABLE III: Energy levels of an infinite (HF)oo chain at the F 
and the X point. Here e{k) represents the Hartree-Fock Bloch 
orbital energies and i?2.2(fc) gives the CO-ADC(2,2) energies 
for the crystal momenta with fc = F, X which were determined 
accounting for an extension of the 2plh and 2hlp configura- 
tions up to third nearest-neighbor cells. All data are given in 
electronvolts. 



trend is observed for (HCl)oo- However, the inclusion of 
nearest-neighbor cells in correlation calculations leads to 
essentially the same displacement of both types of bands, 
i.e., the assumption of a symmetric shift becomes valid 
for long-range contributions as pointed out in Ref. 79. 

Although the quasiparticle bands are shifted appre- 
ciably with respect to the corresponding Hartree-Fock 
bands, they essentially maintain their shape. This ob- 
servation is in contrast to the pronounced reduction of 
bandwidths due to electron correlations which are typi- 
cally found in covalently bonded polymers such as trans- 
polyacetylcnc**" and in covalently bonded crystals such 
as diamond or silicon. ^SiSLS^ 

To investigate the accuracy of the CO-ADC(2,2) 
method, we compare it with theoretical results from 
many-particle perturbation theory which has been ap- 
plied by Liegener and Ladiki^ to the (HF)oo chain cm- 
ploying the DZP basis set. For the top of the valence 
bands and the bottom of the conduction bands at the 
F point, they obtain — 15.10 eV and 3.05 eV, respec- 
tively, in second order of the expansion and — 15.09 eV 
and 3.03 eV in third order. The latter two numbers 
correspond to an increase of the valence band energy 
by 2.48 eV and a lowering of the conduction band en- 
ergy by 0.54 eV with respect to the Hartree-Fock result. 
In contrast, CO-ADC(2,2) yields at the F point an in- 
crease by 2.08 eV and a lowering by 0.33 eV. The over- 
all agreement between CO-ADC(2,2) and many-particle 
perturbation theory is satisfactory. The deviations be- 
tween the data from both methods, as already found for 
the Hartree-Fock band structure, can again be attributed 
predominantly to basis set artifacts and the differences 
in the geometries employed. 

To the best of our knowledge no correlation 
calculations of the quasiparticle band structure of 
(HCl)oo chains have been performed to date. 
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TABLE IV: Energy levels of an infinite (HCl)cx) chain at the 
F and the X point. The labels are chosen as in Tab. lllll This 
time 2plh and 2hlp configurations are being formed by the 
Wannier orbitals in up to second nearest-neighbor cells. All 
data are given in electronvolts. 



VI. CONCLUSION 

This work deals with the quasiparticle band structures 
of infinite (HF)oo and (HCl)co chains, which are deter- 
mined and compared with the energy levels of the con- 
stituting monomers. 

The analysis begins with the examination of the com- 
pounds in Hartree-Fock approximation. The band struc- 
tures of both chains turn out to be quite similar. How- 
ever, the bands of (HCl)oo are shifted towards the band 
gap with respect to the bands of (HF)oo chains. By an- 
alyzing the ionization potentials and electron affinities 
of the monomers, the origin of the individual bands of 
the infinite chains is identified. The higher outer valence 
bands stem from the lone pairs on fluorine and chlorine, 
respectively, and the lower outer valence bands originate 
from the hydrogen bonding. The latter bands also ex- 
hibit less dispersion in (HCl)oo compared with (HF)oo 
due to the weaker bonding. 

Electron correlation calculations are performed to 
determine accurate energy levels of the HF and 
HCl molecules and precise quasiparticle band struc- 
tures of the (HF)oo and (HCl)oo chains. To this end, 
we use the ah initio Green's function theory crystal 
orbital algebraic diagrammatic construction for Wan- 
nier orbitals. The strict second order approximation 
of the off-diagonal Fock matrix elements and the two- 
electron contributions is applied, which is dubbed CO- 
ADC(2,2) methodiS2*2& Setting out from these equations, 
we make the transition to the corresponding method 
for molecules, termed ADC(2,2), by eliminating the lat- 
tice summations. The configuration space of the cor- 
relation calculations is found to converge rapidly; only 
nearest-neighbor cells need to be regarded to determine 
the band gap of the chains within 0.008 eV of the value 
obtained in the largest calculations performed. To assess 
the accuracy of the ADC(2,2) ionization potentials of the 
monomers, we compare with ADC(3) data; a reassuring 
agreement is found. Furthermore, the quasiparticle en- 
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TABLE V: Convergence of the fundamental band gap -Egap 
and the bandwidth of the lower and the upper F 2p va- 
lence band complexes, A.E^2p s'ld ^^F2pi respectively, of 
a (HF)oo chain with respect to the number of unit cells used 
to determine the quasiparticle band structures. "Cells" des- 
ignates the number of unit cells taken into account in CO- 
ADC(2,2) calculations, where zero refers to the original 
Hartree-Fock results. Unity denotes the 2pl/i and 2/ilp con- 
figurations from the origin cell only in correlation calcula- 
tions. Three, five, and seven indicate the additional inclusion 
of nearest-, second nearest- and third nearest-neighbor cells, 
respectively. The top of the valence bands J5top,v and the bot- 
tom of the conduction bands i?bottom,c are both situated at 
the F point. All data are given in electronvolts. 

ergies for the (HF)oo chain compare satisfactorily with 
data from many-particle pertm'bation theory, thus cor- 
roborating the accuracy of the CO-ADC(2,2) scheme. 

The inclusion of correlation effects leads to a pro- 
nounced shift towards the band gap of the quasiparti- 
cle bands with respect to the Hartree-Fock bands. The 
shift of the valence bands is noticeably larger in the 
(HF)oo chain compared with the (HCl)oo chain; the re- 
verse trend is observed for the conduction bands. These 
observations match the behavior of the energy levels of 
the monomers. In contrast to many other compounds, 
e.g., covalently bound polymers^ or crystalS)Si2iiS the 
bandwidth of the chains does not decrease substantially 
due to correlation effects but increases slightly. 

Building on top of this work, one should investigate 
carefully the lower lying, inner valence regime of the 
chains where intriguing many-particle effects are to be 
expectedi^2i^i*Si24 This is indicated by the studies of 
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TABLE VI: Convergence of the fundamental band gap -Bgap 
and the bandwidth of the lower and upper C13p valence 
band complexes, AS^g^ and AS^g^, respectively, of a 
(HCl)oc) chain with respect to the number of unit cells used 
to determine the quasiparticle band structures. The symbols 
are chosen as in Tab. [3 AH data are given in electronvolts. 

oligomers by Deleuze et Q/.j32^33,34^35,36 .^^hich model the 
corresponding infinite periodic chains. They serve to 
investigate the evolution of strong electron correlations 
with an increasing length of the oligomers. Moreover, in- 
termolecular electronic decay processes, e.g., Ref.lssI (and 
references therein) found recently in atomic and molecu- 
lar clusters, including clusters of HF molecules,-^ present 
challenging perspectives. Further emerging effects, when 
making the transition from the monomer of an infinite 
chain, have been discovered in the model oligomer studies 
of Refs. 32 33 34,35.36,87.88,89,90 which deserve a close 
look within an infinite chain treatment. Particularly, the 
examination of these effects in three-dimensional crystals 
offers a fascinating perspective but also represents a great 
challenge. 
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tion matrix, but not the renormalization involvedf^^*^ 
Hence, this nomenclature is imprecise and shall be dis- 
continued. 

In Refs.[ll23, the cc-pVDZ basis set?'^iWi?? was found to 
yield rather inaccurate ground-state binding energies of 
(HF)oo and (HCl)oc chains; they deviated in both cases 
by ~ 0.05 eV from the best value for the binding energy 
[Tab. 4 of R.ef.[lfll|. In this study, however, such accuracy of 
the energy levels would be very good. By these numbers, 
we do not want at all to imply that the computations 
presented here are that accurate but rather indicate the 
significantly lowered demands to the quality of the basis 
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ing Werner and co-workers2^*i- by removing the projected 
atomic orbitals which correspond to the smallest eigenval- 
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Hence, the application of the very fine grained selection 
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